Electron in a tangled chain: multifractality at the small-world critical point 
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We study a simple model of conducting polymers in which a single electron propagates through a 
randomly tangled chain. The model has the geometry of a small-world network, with a small density 
p of crossings of the chain acting as shortcuts for the electron. We use numerical diagonalisation 
and simple analytical arguments to discuss the density of states, inverse participation ratios and 
wave functions. We suggest that there is a critical point at p = and demonstrate finite-size scaling 
of the energy and wave functions at the lower band edge. The wave functions are multifractal. The 
critical exponent of the correlation length is consistent with criticality due to the small-world effect, 
as distinct from the previously discussed, dimensionality-driven Anderson transition. 



Introduction 

The small-world effect occurs in complex networks [l| 
that resemble regular lattices with a small density p of 
"shortcuts" per site 0|. The key observable is the short- 
est path length between two randomly chosen points, I. 
Defining a "thermodynamic limit" by 



oo , p = constant, 



(1) 



where L is the number of sites, one finds that for < 
p « 1 the network has the local structure of a regular 
lattice but some global properties of a random graph Q. 
In fact there is a critical point 0,0] at p = It separates 
regular lattices (p = 0), for which the relative size of 
the shortest path length, l/L, tends to some finite value, 
from lattices with finite density of shortcuts (p > 0) , for 
which l/L — > 0. This quantity shows finite-size scaling 
iflii: for 



i»l,P«l 



(2) 



it depends only on the scaling variable L/£, where £ = 
p~ v is a critical system size diverging algebraically at 
p = 0. A simple renorrnalisation group argument yields 
the critical exponent v [jj]. 

On the other hand it has been accepted for some time 
that conducting polymers must be thought of as complex 
networks 0]. Models have been studied in which a one- 
dimensional Anderson insulator is perturbed by adding 
some additional connectivity in the form of a given den- 
sity of chain "crossings". It is well established that, when 
the density of crossings reaches a finite critical value, 
which depends on the amount of site-energy disorder, 
there is an Anderson metal-insulator transition 



This critical point results from the increase of the ef- 
fective dimensionality of the network. It has also been 
demonstrated [13, E3] in a model having the geometry of 
the Watts-Strogatz (WS) small- world network 0]. 

Interestingly, even models of conducting polymers that 
lack site-energy disorder, so that the state of a single 
wire with no crossings is metallic, display non-trivial be- 
haviour suggestive of proximity to a second-order phase 
transition . But there has been no detailed dis- 

cussion of the nature of this critical point or its location 
on the phase diagram. Here we shall provide evidence 
that it occurs when the density of crossings reaches zero 
and that its origin is the small- world effect. 

The rest of paper is organised as follows. First we 
introduce a simple model that generalises the "winding 
chain" of Ref. [12J and has the geometry of the Newman- 
Watts (NW) network @], with the crossings acting as 
shortcuts. Then we present results of full-spectrum nu- 
merical diagonalisation that demonstrate the small-world 
effect: a qualitative change of the spectrum and wave 
functions on the addition of a small density of shortcuts. 
We then use Lanczos diagonalisation to study in detail, 
for much larger systems, the lower band edge. We find 
finite-size scaling consistent with the critical behaviour 
of the network described in 0, 0| . Simple analytical ar- 
guments are used to characterise the wave functions at 
the critical point. Our conclusions are then laid out. 



Model 

We consider a single particle moving in a tight-binding 
lattice with one orbital per site. The site energy is equal 
on all L sites. Hopping can take place between near- 
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est neighbours or between sites connected by one of TV 
shortcuts. For simplicity we assume the same value of 
the hopping integral, t. The Hamiltonian is 



N 
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where j is a site label and s denotes the shortcut con- 
necting sites i s and j s . We assume periodic boundary 
conditions: \L + 1) = |1). 

Note that t is the only energy scale in the Hamiltonian 
([3]): the system's behaviour is thus completely specified 
by the way in which the TV shortcuts are chosen. One 
can imagine creating these shortcuts by winding an ini- 
tially straight chain, embedded in a higher-dimensional 
space, in a random fashion. Sites that were initially far 
apart may come close so tunnelling of the electron be- 
tween these sites becomes possible. 

The model of conducting polymers in Ref. [12J is ob- 
tained by requiring that such winding is in the form of 
a series of consecutive, non-intersecting loops: i\ < ji < 
*2 < 32 < ■ ■ ■ < ijv < 3N (except for the boundary condi- 
tions, which do not affect our argument). Here we relax 
that restriction: our chain is "tangled" at random so that 
each of 



TV : 



L(L-3) 
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possible shortcuts (excluding only shortcuts joining a site 
to itself or to a nearest neighbour) is present with the 
same probability, 



V 



2p 



(L-3) 



(5) 



The resulting topology is an instance of the NW network 
with critical exponent v = 1 [j?! . 

A similar Hamiltonian, namely the adjacency matrix 
of the WS network, has been considered in [141. Ila|. The 
authors of [HJ solved the time-dependent Schrodinger 
equation, demonstrating the faster spread of initially lo- 
calised wave functions on a small-world network com- 
pared to a regular lattice. In [3| the eigenvalues and 
eigenvectors were obtained, as a means of characterising 
the networks' geometry. Our argument is based on a 
similar calculation. Note, however, that our model has a 
distinct topology, compatible with the simp le p icture of a 
tangled chain — indeed the Hamiltonian of [14 . 1 1 51 ] is not 
of the form ((3j) • Moreover our specific aim is to identify 
critical behaviour near p = [28] . 

In our model the number of shortcuts, N, is a random 
variable following the binomial distribution (Poisson for 
L p). Evidently (N) — Lp, where (...) denotes a 
configurational average. 




Figure 1: (Color online) Density of states (DOS), estimated 
from full-spectrum numerical diagonalisation for (a) R = 32 
realisations, L = 1024 sites and different densities of shortcuts 
per site, p; (b) R = 16, p = 0.1 and different values of L; (c) 
different R with L — 1024 andp = 0.1. All the plots have been 
obtained by drawing histograms using 100 uniformly spaced 
values of the energy from e = —3t to e = +3i. 



Small-world effect 

For a finite size L and a given set of shortcuts (a "re- 
alisation") it is straightforward to diagonalise the Hamil- 
tonian ((3j) numerically [H|- One can then, by averaging 
over a finite number of realisations, R, obtain a coarse- 
grained estimate of the configurationally averaged den- 
sity of states (DOS), p (e) = (J2 U $ ( e ~ £ v)) > f° r a par- 
ticular value of p. Here e v is the energy of the state. 

The main panel of Fig. Q] shows the DOS, estimated 
in this way, for several small values of p. The curve 
for p = 0, given by the well-known formula p(e) = 



e 2 /At 2 j , is also shown. Panels (b) and (c) 
illustrate the dependences of the numerical results on L 
and on R. They suggest that the DOS is well defined 
in the thermodynamic limit ([1]), and that the plot on 
the main panel is representative of it, at this level of 
coarse-graining (higher values of L and R would be re- 
quired to describe the finer detail). At p = 0, the DOS 
has the usual square- root singularities at the band edges. 
For p > 0, the singularities survive in the form of pro- 
nounced maxima in the DOS. This is similar to what was 
described for the model based on the WS network 
Also in common with that model, for p > 1 (not shown), 
the spectral density resembles that obtained by numerical 
diagonalisation of the adjacency matrices of uncorrelated 
random graphs. However the salient features of Fig.[T]are 
two additional peaks outside the conduction band of the 
chain, located at e = ±Vo£. Although this region of the 
spectrum becomes less populated as p — > the position 
of the peaks does not change and they are present for 
any non-zero value of p. They are a qualitative feature 
distinguishing the DOS for p = from that for p > 0. 

To explore further the nature of the states around 
these new peaks we can use numerical diagonalisation 
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Figure 2: (Color online) Energy-dependence of the IPR (left- 
hand axis) and the exponent ct\ characterising its power-law 
dependence on system size L (right-hand axis). The data have 
been obtained by exact numerical diagonalisation for R = 16 
realisations with p = 0.1 shortcuts/site. The IPR vs energy 
curve corresponds to L — 1024 sites. It has been obtained 
by splitting the energy interval from e = — 3t to e = +3t into 
20 sub-intervals and averaging the IPR over each subinterval. 
The 16384 individual values of the IPR in our ensemble are 
also shown. Similar curves have been obtained for different 
values of L, and fits of the averaged IPR to a power-law in 
each energy interval have been used to determine ot\. Two 
examples of such fits are shown in the inset. 



to obtain the inverse participation ratio (IPR; n = 1) 
and higher-order moments of the local density distri- 
bution (n = 2,3,...). These are given by P" 1 (e) = 
Y, jv S(e - e v ) \ip v (j)| 2( " +1) where Vv U) is the wave 

function of the state, evaluated on the j^ 1 site. For 
large enough values of L, we expect the configurationally- 
averaged IPR to be given by a power law: (P^ 1 (e)) = 
L~ ai . The coefficient a\ can be used to distinguish ex- 
tended states (ai = 1) from localised states (evi = 0) and 
states with fractal dimension (0 < a± < 1) [Hj]. Fig. [2] 
shows our results. They suggest that, for p — 0.1, all 
states are fractal, but the wave functions with |e| > 2t 
are much closer to localisation than those with |e| < 2t, 
which are almost extended. Note the strong fluctuations 
of the IPR at particular energies. 

The above results reflect that there are two quite dif- 
ferent types of states. It is easy to show that the real 
eigenfunctions of |(3|) have the form 



fl> (x) = Bie-W'-'d + Cie lK{x ~ Xl) 



(6) 



for xi < x < xi + i, where xi denotes the Z" 1 shortcut ter- 
minal (i.e. xi = i s or j s for some s), ordered by ascending 
site index, and xi + i — xi > 2. The parameter K and the 
corresponding energy €k can only take one of the follow- 
ing forms: K = —in => ex = — 2icoshK; K = k => €k — 
— 2icosfc; K = n — in => eic = 2£cos1ik. Thus in the 
space between shortcut terminals states with |e| < 2t are 
plane waves, with wave number k — cos -1 (— e/2t), while 




Figure 3: Numerically determined single-particle wave func- 
tions for a particular realisation of our tight-binding model, 
with L — 460 sites. Each wave function has been represented 
on a circle, whose radius is its energy, measured from — 3t. 
The innermost wave function in (a) is the ground state. The 
three excited states with energy below —2t have been plotted 
separately for clarity: (b), (c) and (d), in order of increas- 
ing energy. The other wave functions in (a) have energies in 
(— 2t, 0). Finally, the straight lines in the innermost circles 
indicate which sites are joined by "shortcuts". 



those with ±e > 2t decrease or grow exponentially at the 
rate k — cosh -1 (±e/2t). A particular case is the ground 
state of a large chain with a large, single loop, which has 
energy e = — v5t and is exponentially localised at the two 
shortcut terminals, with localisation length k _1 ~ 2 sites 
The case with more than one shortcut terminal is 
more complex. A few wave functions for a particular real- 
isation of our model, with p ~ 0.01 shortcuts per site, are 
shown in Fig. [31 States with e < — 2t are exponentially 
localised on some of the shortcuts. Their energy is very 
close to — VEt, but the degeneracy is broken by localis- 
ing on different sets of shortcuts and by forming bonding 
and antibonding combinations between them. This leads 
to the broadening of the peaks in Fig. [TJ States with 
— 2t < e < resemble plane waves, except that their 
amplitude "jumps" at some of the shortcut terminals. 



Critical behaviour at the lower band edge 

The results presented above indicate a qualitative 
change of behaviour when p becomes non-zero. In par- 
ticular the gradual population of the energy peaks at 
e = ±Vbt is suggestive of a critical point at p = in 
the thermodynamic limit ((TJ) . This further implies the 
universal dependence of all observables on L/£ in the 
finite-size scaling regime |(2]). The latter has been veri- 
fied extensively for geometrical features of the NW and 
WS networks, such as the shortest path length and the 
clustering coefficient Here we are concerned 

with spectral properties and wave functions. To obtain 
data above and below the critical system size £ = p^ 1 
we need to consider very large systems. In this case it is 
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Figure 5: Finite-size scaling of the ground state energy, ob- 
tained as in Fig. [4] The main plot shows the collapse of the 
data for p < 10 -4 on the curve given by Eq. (JS} for values 
of the scaling variable from L/£ = 0.001 to L/£ = 20. The 
two dotted lines indicate eo = — 2t and eo = —v5t. The inset 
shows data for two higher values of p, for comparison. 
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Figure 4: Dependence on the system size L of (a) the energy 
and (b) the IPR [n = 1) and next two higher moments of 
the density distribution (n = 2, 3) at the lower band edge 
obtained by Lanczos diagonalisation and averaged over 2000 
random realisations. The density of shortcuts is p = 5 x 10~ 4 . 
The dotted lines indicate (eo) = — 2t and —VEt (a) and the 
analytical results for p = (b). 



best to use a sparse-matrix numerical algorithm [TJ]. We 
employ which, for real symmetric matrices, is an im- 
plementation of the Lanczos method. For computational 
economy, given that we are now going to examine much 
larger system sizes, we concentrate on the single state at 
the lower band edge. 

Fig. |4] shows, for fixed p, the dependence on L of the 
configurational averages of the ground-state energy eo 
and of the moments of the density distribution for the 
ground state wave function, P ? 7q = IV'o (i)| ■ 
For finite conducting chains at fixed chemical potential, 
the average of eo gives a gross idea of how much the 
energy levels leak out, due to the presence of the short- 
cuts, below the bottom of the original conduction band. 
Our main interest in this quantity, however, is that it 
displays critical behaviour quantitatively consistent with 
the small- world effect (see below). There are four dis- 
tinct regimes: for L <C £ (I) most realisations do not 
contain any crossings, so (eo) ~ — 2t and (P^o) = L~ n ■ 
For L < £ (II) the system is in a crossover region in 
which (eo) decreases and the (P^o) increase with system 
size. For L > £ (III) the energy reaches a plateau and 
stays close to (eo) — y/Et. In this regime most realisa- 
tions have one or more shortcuts, but their terminals are 
quite far apart. The ground-state wave functions look 



like the one in Fig. [3] (a) . Interestingly in this range of 
values of L the (Pno) decrease again. This is because 
although locally the wave function is exponentially lo- 
calised on shortcut terminals, globally the crossings on 
which it chooses to localise can be quite spread over the 
network, and such spread becomes greater as the net- 
work grows in size and more and more crossings become 
available. Finally for L ^> £ (IV) (eo) decreases and the 
(Pno) turn to rise again. This indicates that in the ther- 
modynamic limit (HJ) the lower band edge is at — oo and 
contains localised states, consistent with the Lifshitz tails 
suggested by Fig. [TJ 

Figs. [5] and [6] show (e ) and the (Pno) ( n — 1,2,3) for 
very low values of p, plotted as functions of the scaling 
variable L/£ = Lp. For p < 10~ 4 crossings per site the 
data collapse to a single curve in the ranges of values of 
L/£ shown. The inset to Fig. [5] and Fig. [6] (b) show the 
same data for two higher values of p. Evidently the col- 
lapse is not so good. However there is a range of values 
of Z//£ over which the data depend only on that variable 
which becomes wider as p is lowered. Clearly this widen- 
ing finite-size scaling regime corresponds to the regions 
II and III identified above. 

The defining feature of the finite-size scaling regime is 
that, although the number of shortcuts on the network is 
arbitrary (indeed (N) = £/£), all shortcut terminals are 
well separated. The probability of this is 



sep 



exp [-4p(2Lp- 1) kT 1 ] 



(7) 



(taking N — pL and using that the relevant length scale 
is Ak^ 1 <C L; see the Appendix). Evidently in the ther- 
modynamic limit (Q]) P S ep - * 0. On the other hand the 
finite-size scaling regime is reached when, for a given 
value of the scaling variable Lp, the density of cross- 
ings p is sufficiently low that -Psep ~ 1- In this case 
shortcuts are not close enough to lift the degeneracy of 
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Figure 6: Finite-size scaling of the IPR (n = 1) and next two 
higher moments of the local density distribution (n = 2, 3) 
at the lower band edge, obtained as in Fig. [4] (a) Shows 
the collapse of the data for p < 10 -4 on the curve given by 
Eq. iflQjl for values of the scaling variable from L/£ = 0.1 to 
L/£ = 5. The solid lines correspond to a — 0.48, 0.40 and 0.33 
and fit the n = 1,2,3 data, respectively. The dotted lines 
correspond to a = and a = 1. In (b) data for two higher 
values of p over a wider range of values of L/£ is displayed for 
comparison. 



the ground state and so the only possible energies are 
£o = — 2f (when there are no shortcuts) and eo = —VBt 
(when there are one or more). Since the probability of the 



former is (1 — V) 



(for p -C L), with J\f and V 



defined in Eqs. Q and |(5]) respectively, using VM = pL 
we obtain the following exact analytical expression for 
the finite-size scaling law giving the configurational aver- 
age of the ground-state energy: 

(e > - e- LK {-It) + (l - eT L /«) (-V5t) • (8) 

This is the curve plotted in Fig. \E\ alongside the numerical 
data. Let us now deploy a similar argument for the mo- 
ments (Pnl)- Suppose that, for a given realisation, hav- 
ing N shortcuts, the wave function were exponentially lo- 
calised around a single site xq: ipo (x) ~ e~ K l x-a °l. Then 
the moments at the lower band edge would be given by 



tanh" +1 k/ tanh k (n + 1) . In actual fact the wave 



71,0 



function is localised around 2N' sites, corresponding to 
the terminals of TV' < N shortcuts. Since the sites Xt 
of the terminals are all well separated we can derive (see 
the appendix) 



R 



tanh 



n+l 



71.0 



[tanh/t(rc+ l)2 n iV" 



(9) 



Let us assume that AT' depends only on N and postu- 
late a power law: AT' = N a . Then using the Poisson 
distribution we obtain 



<^o) = E 



(L/0 



N 



tanh 



n+l 



1 



N=l 



tanhK (n + 1) 2 n N 



— (10) 

an v ' 



which gives the finite-size scaling of the moments (P~q)- 
It is important to note however that this expression, un- 
like ([H), contains an adjustable parameter a. Depending 
on its value it can describe very different behaviours. In- 
deed for a — 1 (i.e. requiring that the wave function is 
exponentially localised always on the same fraction of the 
shortcuts) we have 



tanh 



n+l 



2 tanh n (n + 1 



{Ei(L/£)-ln(£/0-7h 



(11) 

where 7 is Euler's constant and Ei (x) is the exponential 
integral. For large £/£, we have (P%o) ~ {L/0 1 (more 
generally (Pno) ~ (L/Q ") corresponding to a state 
which, although exponentially localised at the local level, 
is extended globally over the whole network. Conversly 
for a = (i.e. the wave function is exponentially localised 
always on a single shortcut) we obtain (P^o) ~ constant 
for large L/£, i.e. the wave function is localised in the 
true sense of the word. The dashed lines in Fig. [6] (a) 
represent these two extreme cases. On the other hand the 
data are fitted quite accurately by using three different 
intermediate values of a for n = 1,2,3, respectively. This 
is represented by the solid curves in Fig. [6l It suggests 
that the global structure of the wave function at the lower 
band edge is multifractal in the finite-size scaling regime. 



Conclusion 

In summary we have studied a simple model of con- 
ducting polymers in which chain crossings act as short- 
cuts for the electrons, so they move in a small-world 
topology. Our results are based on numerical diagonal- 
isation (full spectrum for systems of size L ~ 10 3 and 
Lanczos for L <~ 10 5 ) and some simple analytical argu- 
ments. We have seen that a small density of shortcuts p 
leads to qualitative, but continuous changes in the DOS 
and wave functions. The former take the form of new 
peaks appearing at well defined positions e = ±y/Et out- 
side the initial conduction band, |e| < 2i. At these peaks, 
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the wave functions have a non-trivial structure. We have 
investigated in detail the lower band edge, confirming 
through finite-size scaling that there is a critical point at 
p = 0. The critical exponent of the correlation length is 
consistent with criticality due to the small- world effect. 
We have derived analytical expressions for the scaling 
laws of the energy and the IPR and higher moments of 
the density distribution. The latter expressions contain 
an adjustable parameter a describing the global structure 
of the wave function, which we find to be multifractal. 

We end by making two additional remarks. Firstly, 
the critical point described here is quite distinct from the 
dimensionality-driven Anderson transition considered be- 
fore 0, 0, 0, 11 1 . The latter occurs at a finite crit- 
ical density of crossings p c (W) > 0, which depends on 
the amount of site-energy disorder W. Since real highly 
conducting polymers are intrinsically disordered the 
small-world critical point at p = is not directly acces- 
sible to experiments. Nevertheless it may make a sig- 
nificant contribution to conductance fluctuations in the 
metallic state. Secondly, the multifractality of the wave 
functions suggests that the adjacency matrices of small- 
world networks might be describable by the power law 
banded matrix ensemble of random matrix theory. 
This intriguing possibility remains largely unexplored 

We thank P.R.A. Campos, K. Capelle, E.C. Carter, 
J.M.F. Gunn, C. Hooley, V.E. Kravtsov, I.V. Lerner, 
B.N. Narozhny, L.N. Oliveira, A.J. Schofield, R.A. Smith 
and I.V. Yurkevich for useful discussions. JQ thanks the 
Abdus Salam ICTP for hospitality during the prepara- 
tion of part of this manuscript. We acknowledge financial 
support from CAPES (VLC), FAPESP and The Lever- 
hulme Trust (JQ). In addition, JQ gratefully acknowl- 
edges an Atlas fellowship granted by CCLRC in associ- 
ation with St. Catherine's College in the Univeristy of 
Oxford. 



APPENDIX 

In this Appendix we provide some simple derivations 
omitted in the main text. Firstly, let us prove that the 
probability P se p of all the 2N shortcut terminals being 
well separated from each other is given by Eq. ([7]). Each 
terminal is found at a random location along the chain. 
Let Xt, with t — 1,2, ...,2N, be the site indices of those 
locations. As noted in the main text, the ground-state 
wave function decays exponentially as e - K \ x ~ x t\ i n the 
space between terminals. After selecting the location of 
the first terminal, the probability for at least one of the 
2N — 1 remaining terminals to be within 2k _1 of x\ (the 
factor of two is arbitrary and does not affect the argu- 
ment) is (4k- 1 /L)(2N - 1). Thus the probability that 
the first terminal is separated form all the other termi- 
nals (in the sense that it is within two decay lengths of 



any one of them) is 



P 1 = l--^-(2JV-l). 



(12) 



Now consider the second terminal, located at X2- The 
probability of the remaining 2N — 2 terminals being well 
separated from this one, in the same sense as above, is 
1 - {4k- 1 /L)(2N - 2). So the probability of both the 
first and second terminals being well separated from all 
the others is 



P 2 = 



4k~ 1 



1 - ^— (2N - 2) 



(13) 



At this point we can already see that the probability of 
all the terminals being well separated from each other is 
given by 



2JV-1 r 

-Psep = P2N-1 = Yl 
whence, taking logarithms, 



1_ (2N-j) 



(14) 



log (Psep) 



4k~ 



2JV-1 



E 2N - 



4k- 1 N 



3=1 



(2N- 1). 
(15) 

Considering that the averaged number of shortcuts is pL, 
Eq.® follows. 

Let us now prove Eq. J9]) for the moments of the density 
distribution of the ground-state wave function. With N 
well separated shortcuts, let us suppose the ground-state 
wave function has peaks at 2N' shortcut terminals (N' < 
N) , decaying exponentially from each maximum again as 

e -K\x-X t \^ rp^g momen t p-l Jg gj ven D y 



Pnl = E^0(i)| 2 " 
3=1 

2N' 

w \A\ 2n+2 J2 



\A\ 



2n+2 



J- ( e -«|j|)2n+2 
j=-oo 

2N' 



tanh [{n + 1) k] 



(16) 



where A is a normalization constant. The latter can be 
found from Pq 1 = 1. It is 



2 _ tanh(/c) 
2AT' ' 

which substituted in Eq. lfT6|l yields Eq. ([9]). 



(17) 
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